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ABSTRACT 

Massive planetary cores (~ 10 Earth masses) trigger rapid gas accretion to form gas giant 
planets such as Jupiter and Saturn. We investigate the core growth and the possibilities for cores 
to reach such a critical core mass. At the late stage, planetary cores grow through collisions with 
small planetesimals. Collisional fragmentation of planetesimals, which is induced by gravitational 
interaction with planetary cores, reduces the amount of planetesimals surrounding them, and 
thus the final core masses. Starting from small planetesimals that the fragmentation rapidly 
removes, less massive cores are formed. However, planetary cores acquire atmospheres that 
enlarge their collisional cross section before rapid gas accretion. Once planetary cores exceed 
about Mars mass, atmospheres significantly accelerate the growth of cores. We show that, taking 
into account the effects of fragmentation and atmosphere, initially large planetesimals enable 
formation of sufficiently massive cores. On the other hand, because the growth of cores is slow 
for large planetesimals, a massive disk is necessary for cores to grow enough within a disk lifetime. 
If the disk with lOOkm-sized initial planetesimals is 10 times as massive as the minimum mass 
solar nebula, planetary cores can exceed 10 Earth masses in the Jovian planet region (> 5 AU). 

Subject headings: planet and satellites:formation 



1. Introduction 

Gas giant planets such as Jupiter and Sat- 
urn form in gaseous disks. In the core-accretion 
model, the accretion of planetesimals produces 
cores of giant planets. Once a core reaches a 
critical mass ~ 10 Earth masses, it can rapidly 



accre t e gas to form a gas giant pl anet (Mizuno 



1980; iBodenheimer fc Pollack! Il986t llkoma et al 
2000). Gas giants must form within the lifetime 
of gaseous disks (< lOMyr). 

For km-sized or larger planetesimals, gravita- 
tional focusing enhances their collisional cross sec- 
tions, resulting in a high collision probability for 
low relative velocities. Relative velocities of large 
bodies are kept lower than those of small ones due 
to dynamical friction. A combination of gravita- 



tional focusing and dynamical friction brings rapid 
growth of lar ge bodies, which is referred to as run- 
away growth (jWetherill fc StewartHl989|) . Eventu- 
ally, the runaway growth generates a small popu- 
lation of large bodies called planetary embryos. 
Planetary embryos keep their orbital separations 
and hence grow through collisions with surround- 
ing remnant pl anetesimals more slow ly than in the 
runaway mode (jKokubo fc I da 1998). This regime 
is called oligarchic growth. 

iKobavashi etaH (|2010h pointed out that the 
oligarchic growth halts due to fragmentation of 
planetesimals. In the oligarchic growth, the rel- 
ative velocities of planetesimals are controlled by 
the viscous stirring of embryos and gas drag. As 
embryos grow, the velocities of remnant planetes- 
imals are increased so greatly that collisions be- 
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tween planetesimals become destructive. Such 
collisions eject numerous fragments, which col- 
lide with each other to produce further smaller 
bodies. Planetesimals are therefore ground down 
through such successive collisions (collision cas- 
cade). The random velocities of small bodies are 
strongly damped by gas drag and thereby the colli- 
sional cascade no longer occurs for fragments with 
radii < 1-10 m. In the end such fragments drift 
inward due to gas drag and are lost around em- 
bryos. The collisional cascade combined with the 
loss of fragments reduces the soli d surface density 



Kobavashi et al 



and h ence final embryo masses 
(2010) showed that the final embryo masses are as 
small as Mars mass for 1-100 km-sized planetesi- 
mals in t he minimum m ass solar nebula (hereafter 



MMSN; lHavashil 119811 ). Large planetesimals, 



which are relatively hard to be broken collision- 
ally, and a massive disk produce massive final em- 
bryos. However, collisional fragmentation makes 
it difficult to form giant planets along the lines 
of the core-accretion model; starting from 100 km- 
sized planetesimals, planetary embryos can reach 
the critical core mass for gas accretion only inside 
3-4 AU in a disk that is 10 times more massive 
than the MMSN model. 

The motion of fragments < 1 m is coupled with 
gas. The drif t timescale of such fragmen ts are rel- 
atively long. iKenvon fc Bromley ( 20091 ) proposed 
that embryos may accrete a large amount of such 
fragments. However, the strong gas drag in the 
Stokes regime is dominant for fragments < 100 m 
and damps the relative velocities to halt collision 
cascade at 1-10 m as mentioned above. Therefore, 
only a small amount of coupled bodies are pro- 
duced a nd hence they hardly co ntribute to embryo 
growth ( Kobavashi et al. 2010h . 

Many authors have investigated embryo growth 
with i V-body, statistical, a nd h y brid si mula- 
tions dKokubo fc Ida! 119961 Il998l l2000l [2 002: 



Inaba et al.lll999ll200lll2003t[Weidenschilling et al 



1997t Weiderischillinel2005ll2008l^ jKenvon fc Bromley 
2004j2008l:IChambersl200dl2008l:lKo"bavashi et al " 



2010). Although providing most accurate dynam- 
ical results, TV-body simulations have difficulty 
in producing numerous fragments and following 
their fate. The fragmentation effect on embryo 
growth has thus not been treated in detail in spite 
of its importance. Recently, Levison et al. (2010) 
included fragment production in their TV-body 



simulation. However, it is still difficult to treat 
fragment-fragment collisions. Such successive col- 
lisions are essential in the collision cascade (e.g., 
Kobavashi fc Tanaka 2010l ). Therefore, statisti- 



cal simulations are a better method to accurately 
investigate planet formation with fragmentation. 

In the statistical simulation, the collisional 
mass evolution of bodies is calculated within a 
"particle-in-a-box" approximation. Bodies have 
horizontal and vertical components of random ve- 
locity relative to a circular orbit that are deter- 
mined by their eccentricities and inclinations, re- 
spectively. These velocities are changed by gravi- 
tational interactions between the bodies and hence 
affected by their mass spectrum, while the colli- 
sion rates between the bodies depend on the ve- 
locities. Therefore, the coupl ed mass and velocity 

Steward 



evolu tion needs to be solved (jWetherill 
While 
its 



1993) 



the statistical method has ad- 
weak point is the inability to 



vantages, 

track the individual positions of planetesimals. 
However, progress in p l anetar y dynamic theory 
dGreenzweig fc Lissauer! Il992l: llda fc Nakazawa 
1989[l0htsukill999l : IStewart fc Idal2000HOhtsukret all 
2002h has h elped to overcome th is pro blem. 
For example. iGreenzweig fc Lissauerl (|1992) and 
Ida fc Nakazawal (|1989l ) provided detailed ex- 



pressions for the probability of collisions be- 
tween planetes i mals o rbitin g a central sta r, while 
IStewart fc Idal (|2000l ) and lOhtsuki et all (|2002h 



derived improved equations for calculating the 
evolution of random planetesimal velocities caused 
by gravitational interactions. Finally, it has 
been shown that the recently developed statis- 
tical codes can describe some aspects of the plan- 
etary accumulation processes with the same ac- 



curacy as iV-body sim ulations ([Inaba et al.l 12001 
Kobavashi et al.H2010l ). 



Since the timescale of collision cascade strongly 



affect s the final mass of planetary embryos (jKobavashi et al 
201dh . fragmentation outcome models are es- 
sential for embryo growth. Collisional frag- 
mentati on includes several uncer tain parame- 
Kobavashi fc Tanakal (|2010l ) constructed 



ters. 



a simple fragmentation model w hich is consis- 
tent w ith laboratory expe r iments dFuiiwara et al 



1 19771 : iTakagi et al.l 1 19841: iHolsappld Il993l) and 



hydro dynamical simulations ( Benz fc Asphaud 
1999) and analytically clarified which parame- 
ters are essential. They found that the mass de- 
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pletion due to collision cascades is sensitive to 
the total ejecta mass yielded by a single colli- 
sion, while it is almost independent of the mass 
of the largest ejecta fragment and the size dis- 
tribution of ejecta over a realistic parameter re- 
gion. Furthermore, fragmenting collisions are sub- 
divided into two types, catastrophic disruption 
and cratering (erosive collision). Although some 
studies neglected or underestimated the effect of 



Il969t 

1994tlWetherill fc Stewartll 19931 ; llnaba et al.lbool 
Bottke et al. 2005t h Kobavashi & Tanaka showed 
that cratering collisions make a dominant contri- 
bution to the collision cascade. 

A planetary embryo larger than ~ 10~ 2 Earth 
masses acquires a tenuous atmosphere of gas from 
the disk. Fragments are captured by the atmo- 
sphere even if they do not collide directly with 
the embryo, implying that the collisional cross 
section of the embryo is enhanced. This ef- 
fect adva nces the growth of M ars-mass or larger 
embryos (jlnaba fc Ikomal [2003) . Embryos with 
the atmospheres accrete fragments prior to their 
drift inward and can acquire more than 10 Earth 
masses starting fr om lOkm-sized planetesimals 



cratering (IDohnanvil Il969t IWilliams fe Wcthcrill 



( Inaba et al.l 12003 ). However, erosive collisions 



and initial plan etesimal sizes strongly affect final 



embryo masses (|Kobavashi et al.l f2010) 



This paper investigates the embryo growth tak- 
ing into account erosive collisions and embryo's 
atmosphere. Although growing embryos may fall 
into a central star d ue to the type I migration (e.g., 
Tanaka et al.ll2002h . we neglect the migration here. 
We perform both analytical studies and statistic 
simulations, which extend those of Kobayashi et 
al. (2010) by including atmospheric enhancement 
of embryo growth. The goal is to find out what de- 
termines embryo growth and whether an embryo 
can reach the critical core mass. We introduce 
the theoretical model in Section [2] and derive final 
embryo masses taking into account atmosphere in 
Section [3] In Section [4j we check solutions for fi- 
nal masses against the statistical simulations. Sec- 
tions [5] and [6] contain a discussion and a summary 
of our findings. 



2. THEORETICAL MODEL 
2.1. Disk Model 

We introduce a power-law disk model for the 
initial surface mass density of solids E Sj o and gas 
So- o such that 



Sgas,0 



/iccEl (lALl) ' g< 
/gasSl (lAlj) g 



(1) 
(2) 



where a is a distance from a central star, Ei is 
the reference surface density at 1 AU, and q is the 
power-law index of the radial distribution. The 
gas-dust ratio / gas = 240 (Hayashi 1981). The fac- 
tor /i C e that represents the increase of solid density 
by ice condensation beyond the snow line a lcc is 
given by / icc = 1 (a < a icc ) and 4.2 (a > a icc ). In 
the MMSN model, Si = 7.1 gem" 2 and q = 3/2. 
If the disk is optically thin, 



2.7 



Li 



1/2 



AU, 



(3) 



where L» and L© are the luminosities of the cen- 
tral star and the sun, respectively. In reality, disks 
may be optically thick even after planetesimal for- 
mation. However, we assume Equation (j3|) for sim- 
plicity. 

2.2. Fragmentation Outcome Model 

Kobavashi fc Tanakal (|2010l ) showed erosive col- 
lisions to dominate the collision cascade. We 
should take into account such collisions prop- 
erly. We assume that fragmentation outcomes are 
scaled by the impact energy and hence the total 
ejecta mass m e produced by a single collision be- 
tween mi and ni2 is given by a function of the di- 
mensionless impact energy <f> — m\m2V 1 j2(m\ + 
to 2) 2 <3d> wnere v is the collisional velocity be- 
tween mi and mi and is the specific en- 
ergy needed for m fi = (mi + mp)l2. Following 
Kobayashi fc Tanakal ( 20101 ) and iKobavashi et al 
J20irj|), we model 



mi + ni2 



(4) 



Inaba et al.1 (|2003h derive d m R from the frag- 



ment model developed by IWetherill fc Steward 
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1993 ) with a value of QJj found bv lBenz fc Asphauel 



1999) for ice. Fig.Q]shows t heir model and Equa- 
tion (|H) . As discussed in iKobavashi fc Tanakal 
( 20101 ). most of the laboratory experiments and 



the hydrodynamic numerical simulations of colli- 
sional disruption showed m P , not to have a discon- 
tinuity at (ft = 1 (IH ouscn e t al.lll991t iTakagi et al 



19841 : iBenz fc Asphaud" Il999t ). Therefore, Equa 
tion (j4j) includes erosive collisions ((ft < 1) more 
accurately. 



+ 



0.001 




0.01 



Fig. 1. — The total ejecta mass m e produced by 
a single collision with mi and 7712 , as a function 
of the dimensionless energy (ft = mim-iv 2 /2(m,i + 
ni2) 2 Q^,- The solid line indicates Equation (j4]). 
For reference, the dotted lines are sh own for the 
fragment model of llnaba et al.l (]2003f ) with mi = 
10 3 m 2 = 4.2 x 10 20 g. 



The critical energy QJj is given by 

f r V B ( r V s 
Qd = Qos \ 7777 ) + QogPp - I + C K 



2Gm 



1 cm 



1 cm 



where r and m are the radius and mass of a body, 
p p is its density, and G is the gravitational con- 
stant. The first term on the right-hand side of 
Equation ([5|) is dominant for r < 10 4 -10 5 cm, 
the second term describes QJj of r < 10 7 cm, 
and the third term contr o ls Qp for the larger 



bodies. IBenz fe Asphauel (|1999n performed the 



hydrodynamical simulations of collisional disper- 
sion for r = 1-10 7 cm and provided the values of 
Qos, /3 S , Qog, and /3 g . For r > 10 7 cm, is purely 



determined by the gravitational binding energy, 
being independent of material properties. The 
collisional simul ation for gravitational aggr egates 
yields C gg ~ 10 (|Stewart fc Leinhardtl [20091) . 



2.3. Enhancement Radius by Atmosphere 

Once a planetary embryo has grown larger than 
the Moon, it acquires an atmosphere. It helps the 
accretion of planetesimals or fragments onto an 
embryo; small bodies are captured by the atmo- 
sphere of the embryo. 

Inaba fc Ikomal ( 2003 ) provided an analytical 



model for a density profile of the atmosphere. We 
consider the atmosphere at a distance R c from 
an embryo center. We assume that R c is much 
smaller than that at the outer boundary of the at- 
mosphere and that its temperature is much higher 
than that at the boundary. The at mospheric den- 
sity <9 a , is then p roportional to R~ 3 ( Mizunol 1980; 
Stevenso'nl 19821 ). Applying the temperature T n eb, 
pressure P n cb, and density p nc b of the nebula in the 
disk midplane as those at the outer boundary of 
atmosphere, the density profile of the atmosphere 
around an embryo with mass M is given by 



Pa {R e ) _ 16Tra SB GMT^ h fGMp neb 



Pncb 



3ftZ/ -P nc b 



, 4-Pneb-^e 



(0) 



where k is the opacity of the atmosphere and ctsb 
is the Stephan-Boltzmann constant. The plane- 
tary luminosity L e mainly comes from the accre- 
tion of bodies. We approximate 



GM dM 



(7) 



where R is the embryo radius. To validate the as- 
sumption of p a oc R~ 3 , we will ap ply the complete 



W model by Inaba fc Ikomal ( 2003 ) to our statistic 
simulation and compare our analytical solutions 
with the statistical simulations in Section 0] 

When a body passes by a planetary embryo 
with an atmosphere, the embryo can accrete the 
body without direct collision due to the atmo- 
sphere. The relative velocity between the body 
and the embryo at infinity is determined by the 
eccentricity e of the small body; it is given by ev\ 
with the Keplerian velocity = y GM*/a and 
M* being the mass of a central star. The rela- 
tive velocity is typically smaller than the surface 
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escape velocity of the embryo during the embryo 
growth. If the orbital energy of a body is suffi- 
ciently reduced by the atmospheric gas drag, the 
body is captured by the embryo. The maximum 
rad ius r of bodies capture d at distance R e is given 
by (jlnaba fc Ikomall2003l) 



9ah M Pa,(Re 



6 



(8) 



where h M = (M/3M,) 1 / 3 is the reduced Hill 
radius of the embryo and e = ejhu- Equa- 
tion (|p is derived under the two- b ody approxima- 
tion. iTanigawa fe Ohtsuki (120101) confirmed that 



Equation © is valid in the case where the three- 
body effects are included. 

Equation ([8]) means that R e is the effective col- 
lisional radius of an embryo for bodies with radius 
r. The enhanced radius of the embryo with atmo- 
sphere is thus derived from Eqs. (jl])-© as 



R c 
~R 



FM 8 / 9 

TO 1 /9Ml/3 : 



where 



F 



^ 2 aas B T n 4 eb/ 4 b G3 



-,1/3 



(e 2 +6)(3M»)V3 K p4 ob 



(9) 



(10) 



The enhancement factor R e /R given by Equa- 
tion © is shown in Fig. [3J where the power- 
law density profile given by Equation (J6j) is 
compared with a mo re realistic profile given by 
Inaba fc Ikomal (|2003h . As we discuss later, plan- 
etary embryos mainly grow through collisions with 
planetesimals of the initial size or with fragments 
of radius r ~ 10 m. The enhancement factor cal- 
culated with Equation (JSJ reproduces well the 
more realistic one for km-sized or larger planetesi- 
mals, but Equation (0) significantly overestimates 
R c /R for fragments. However, since the accre- 
tion rate due to collision with such fragments has 
a weak dependence on the enhancement factor 
(oc (Rc/R) 1 / 2 ; see Equation (|30| ). this discrep- 
ancy produces insignificant errors. 

3. FINAL EMBRYO MASS 
3.1. Isolation Mass 

Planetary embryos can grow until they have 
accreted all planetesimals within their feeding 



ft 




Fig. 2. — The ratio of the enhanced radius of plan- 
etary embryo to its physical radius with M = Mq , 
p p = lgcm~ 3 , and M = 1 x 10 -6 M e /yr for 
e = 4 in the MMSN disk around the star with 
mass Mp, . The ratios are calc ulated by the formu- 
lae of llnaba fc Ikomal (|2003l ) for the opacity ob- 
tained from Equation ([31?]) with the grain deple- 
tion factor / = 10~ 4 1 (solid lines) and by Equa- 
tion for the constant opacity n = 0.01 cm 2 g _1 
and lcm 2 g _1 (dotted lines). 



zones. The width of a feeding zone is given 
by the orbital separation of neighboring em- 
bryos, 6(2M/3M») 1 / 3 a, where b ~ 10 is the 
separation measured in th eir mutual Hill radii 
(jKokubo fc Idal 1200(1 1 2002). The maximum mass 
or "isolation mass" is M iso = 27ra 2 (2A/ iso /3Af») 1 / 3 6S Si0 . 
It can be expressed as 



M iso = 




where is the Earth mass and M Q is the 

solar mass. The planetary embryo mass ap- 
proaches the_js^lation mass if fra gmentation is 
ignored (|Kokubo fc Idal 1200(1 EEj ). However, if 
fragmentation is included, the embryo mass can 
reach only about Mar s mass for a MMSN disk 
(|Kobavashi et alj|2010f ). 
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3.2. Planetesimal Accretion 



As shown by iKobavashi etaH (|2010h . a plan- 
etary embryo accretes planetesimals with masses 
comparable to original ones or fragments resulting 
from collisional grinding of planetesimals. In the 
former final embryo mass is determined by 

the equilibrium between the accretion of planetes- 
imals and their removal due to collisional grind- 
ing. In the latter case, an embryo can grow un- 
til frag ments are deplet e d by the gas drag. Fol- 
lowing IKobavashi et al. ( 2010l ) , we want here to 
derive final masses determined by the accretion 
of planetesimals in the case with atmospheric en- 
hancement, while we treat the fragment accretion 
in Section [B~3l 

At the oligarchic stage, embryos mainly grow 
through collisions with planetesimals that domi- 
nate the surface density. The growth rate of an 
embryo with mass M is given by 



dM 



C a 



~'E,a 2 h 2 



Pr.nl ) 



M V* col 



(12) 



where JIk is the Keplerian frequency and C acc is 
the correction factor on the order of unity. The 
dimensionless collision rate (P co i) is formulated as 
a function of the eccentricities e and inclinations 
i of bodies accreted onto the embryo. We assume 
e = 2i in this analysis. 

Embryos have a consta nt ratio of their sepa - 
rations to their Hill radii ( Kokubo fc Ida 1998). 
When the ratio decreases as embryos grow, rela- 
tively smaller embryos are culled and thereby re- 
maining embryos keep the ratio constant. Sup- 
posing the cull occurs instantaneously, the growth 
rate of embryos due to the cull is estimated to be a 
half of that fro m planetesimals. Therefore, we set 
Cacc = 1-5 (e.g.. IChambersll200d ; IKobavashi et al 



201(1 ). 



For kilometer-sized or larger planetesimals, 
their eccentricities e are controlled by the em- 
bryo stirring and gas drag. The stirring rate 
is written as de 2 /dt — um d 2 h\j(Pys)^k, where 
riM is the surface number density of embryos and 
the dimensionless stirring rate (Pvs) is given by 
(P vs ) = C vs h 2 M ln(A 2 + l)l e 2 with C vs = 40 an d 
A = 5e 3 /96 for e » h M (jOhtsuki et all 120021) . 
Although ln(A 2 + 1) in (Pvs) is weakly dependent 
on e, we adopt, in this analysis, In (A 2 + 1) ~ 8, 
with which va l ue we can reproduce the formula of 
Ohtsuki et ail (|2002h for e = 3-10. The gas-drag 



time r is characterized as ( Adachi et al. 19761) 

2m 



-rrr 2 C D p neh v k 



(13) 



where the dimensionless gas drag coefficient Cd = 
0.5 for km-sized or larger planetesimals. It should 
be noted that r is the stopping time due to gas 
drag only when the relative velocity u between 
gas and a body is equal to the Keplerian ve- 
locity; hence r is almost always much shorter 
than the stopping time for realistic relative ve- 
locities. The e-damping rate due to gas drag is 
given by de 2 Jdt = — C gas e 3 /r with C gas = 2.1 
(|lnaba et al.l 120011 ). Using um — (2iraSa)~ 1 with 
the orbital sepa ration of neighboring embryos of 
6a = 2 1 ^h M ab (|Kokubo fc Idal l2000h and equat- 



ing the stirring and damping rates result in the 
equilibrium eccentricity: [J 



C VS m(A 2 + l)n kT 



2 4 / 3 7T&C*, 



gas 



1/5 



(14) 



Since we roughly estimate e ~ (rf^) 1 ' 5 from 
Equation (j!4|) . the eccentricities of the kilometer- 
sized and larger bodies are larger than Km, ac- 
cording to the assumption e 3> Km- 

Taking into account the enhancement due 
to the atmosphere, the dimensionless colli- 
sional probability for e Km is given by 



dGreenzweig fc Lissauerl Il992t llnaba et al 
Inaba fc Ikomall2003l) 



2001: 



(Pool) 



C co \R R c 
P ~R 



(15) 



where C col = 36 and R = R/ah M = (9M 4 ,/47r / 9 p ) 1 / 3 /a. 
Inserting Eqs. © and (fT5)) to Equation ((T2)) . we 
obtain M as 



^ = A ca M 7 ^/\ 
dt 



where 



C acc a 2 C co iRFQ k 
(3Af„) 2 / 3 g 2 m 1 /9 



3/4 



(16) 



(17) 



J |lda fc Makind dl993T) and iThommes et al.l l|2003fl pre- 
sented a similar equation from the stirring timescale de- 
rived by llda fc Makind Jl993f) ■ We apply the formula of 

lOhtsuki et al .1 (20021 ). which weakly depends on e through 
ln(A 2 + 1). However, since we adopt a constant value for 
ln(A 2 + l) in this analysis, there is no substantial difference 
between their and our treatment, except for the definition 
of the coefficient for the viscous stirring. 



6 



As embryos grow, destructive collisions be- 
tween planetesimals are induced by the stirring 
of embryos and generate a lot of small fragments, 
which produce further small bodies through mu- 
tual collisions. Since very small bodies re- 
sulting from successive collisions are rapidly re- 
moved by the gas drag, the collision cascade re- 
duces the surface density of solids. In the colli- 
sion cascade, collisional fragmentation dominates 
the mass flux along the mass coordinate. Since 
the mass flux is independent of mass in a steady 
state, the mass distribution of fragments follows a 
power law and the power-law exponent a is given 
by ot = (11 + 3jp)/(6 + 3jp ) for e 2 /Qfj oc mr p 
( Kobavashi fe Tanaka 2010l ). The steady-state 
mas s flux determines the surface density reduction 
as (IKobavashi fe Tanaka 20101 Kobavashi et al 
2010h 



dt 



where 



Sl23(«) 



-B ca E 2 M 2(Q - 1)/3 . 
(2-a) 2 O k s 123 (a) 

TO 1 / 3 
~2 2 

2(3Af*) 2/3 gj 



>ln ■ 



2-6 r (1 



(18) 



(19) 



(20) 



2/3 - 

and ho = l.lp p 1 . For the derivation of Equa- 
tion (fTB"]) . we apply the fragme nt ation outcome 



model of IKobavashi fe Tanakal ( 2010l ): ejecta 
yielded by a single collision between m\ and mi 
are characterised by their total mass m G and their 
power-law mass spectrum with an exponent 6 be- 
low the mass — e[m\ + m-i)§l{\ + </>) 2 , where 
e < 1 is a constant. The E s reduction rate is in- 



sensitive to e and 6 ( Kobavashi fe Tanaka 2010). 
We set 6 = 5/3 and e = 0.2 in this paper. 

Dividing Equation (|16|) by Equation (|18|) and 
integrating, we obtain the relation between the 
embryo mass M and the surface density E s : 



G 



4a 



M (4a-5)/6 _ M( (4a-5)/6 



= 4(s s - i / 4 _ Es: y 4) ^ j (21) 



where AIq is the initial embryo mass. Note that 
the derivation of Equation (f2~T|) assumed that the 
planetesimal density reduction is caused by col- 
lisional grinding, but the planetesimal accretion 
onto embryos significantly contributes to the E s - 
reduction when the surface density of planetesi- 
mals, E s , is much smaller than that of embryos, 
Mum- When an embryo reaches a final mass 
M ca , E s may be described as C-z a M ca riM with a 
constant Cs <C 1; hence 



E s 
E s ,o 



{M c: 



2/3 



(22) 



For C Eb ~ 0.1, a final mass is co nsistent with sim- 
ulations ( Kobavashi et al. 2010l ). We thus set 
Cs 3 = 0.1 to derive a final mass. From Eqs. (|2T|) 
and (|22p . we obtain a final embryo mass 



2(4a-5)A ca C- 1/4 S- 1 / 4 M i 1 s / 6 

3i?ca 



3/2(a-l) 



(23) 



Here, we assume M ca » M . 

For kilometer-sized or larger planetesimals, 
Qr> = < 9ogPp r ' 3g with constants Qog and (3 g . We 
apply Qty = 2.1 erg cm 3 g~ 2 and (3 g = 1.19 for ice 
( Benz fe Asphauglll999h and e 2 3> 6, and Equa- 
tion (|2"3"|) can then be re- written as 



M ra = 1.8 x 10~ 



5AU 



2.8 



4 x 10 20 g 



0.63 



2.1 erg cm 3 g~ 2 

/gas El 

1.7 x 10 3 gem' 



1.5 



0.01 gem" 



-0.51 



(24) 



Since planetesimals grow before planetesimals' 
fragmentation starts, planetesimal mass m is 
slightly larger than initial planetesimal mass mo- 
Kobayashi et al. (2010) showed that planetesimals 
mainly accreting onto embryos have m = IOOtoo- 
For mo > 10 23 g (tq > 3 x 10 3 km), final embryo 
masses exceed 10 M e at 5 AU in a MMSN disk, 
but embryos cannot reach it within a disk life- 
time due to their slow growth. The final mass 
M ca is independent of E Sj o, while high E g ,o in- 
creases M ca because gas drag highly damps e. For 
Ei = 71gcnT 2 (lOxMMSN), initial planetesi- 
mals with r > 50 km can produce an embryo 
with 10 M ffi at 5 AU. 
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For comparison, we also show the final mass 
M c in the same situation but neglecting the atmo- 
sphere (Kobayashi et al. 2010): 



0.63 



M C = 0.10 

'l n(S s ,o/S s ) 
4.5 



in 



AU/ V 4xl ° 20 ! 

1.21 



0.48 



O.S!) 



1.7 x 10 3 gem" 



2.1 erg cm 3 g~ 2 / 

1.21 

M e , (25) 



where ln(E Sj o/S s ) — 4.5 is estimated from Equa- 
tiondU]) with C Ss =0.1 for M = 0.1M @ in the 
MMSN model. The collisional enhancement due 
to the atmosphere is inefficient for m = 4 x 10 20 g; 
M ca < M a . If m > 4 x 10 22 g, the atmosphere 
contributes to embryo growth. 

3.3. Fragment Accretion 

As described above, planetesimals are ground 
down by collision cascade and resulting small frag- 
ments spiral into the central star by gas drag. 
In the steady state of collision cascade, the sur- 
face density of planetesimals is much larger than 
that of fragments. However, when the grinding of 
planetesimals is much quicker than the removal of 
small fragments by gas drag, fragments accumu- 
late at the low-mass end of collision cascade and 
determine the total mass of bodies. Embryos then 
grow through the accretion of such fragments. 

The specific impact energy between equal-sized 
bodies, e 2 w 2 /8, should be much smaller than Qjj 
at the low-mass end; thus the typical fragments at 
the low-mass end have 



2 2 

e v k 



(26) 



where Cl ~ 1 is a constant. Although lKobavashi et al 
(|2010h used Cl = 1 to determine the typi- 
cal fragment mass, we apply Cl = 0.5 to cor- 
rect a mistake of factor 2 in their e 2 . Such 
small fragments feel strong gas drag in Stokes 
regime; Cd = 5.5cl g /ur, where c is the sound 
velocity and Z g = lg,o/p g is the mean free path 
of gas molecules w ith l gy0 — 1.7 x 10~ 9 gcm~ 2 



(Ad achi et al. The eccentricities of frag- 



ments at the low-mass end are much smaller than 
h,M and T), where r/ — (v^ — w gas )/wk is the de- 
viation of the gas rotation velocity v gas from the 



Keplerian velocity. The dimensionless viscous stir- 
ring rate i s given by {jVs} = (-Pvsjow) = 73 for 
e <C Hm (jOhtsuki et al.l 120021 ) and the damping 
rate is expressed as de 2 /dt = —2rje 2 /T for e <C r] 
(Ad achi et al. I ll976l) . The equilibrium eccentricity 



between stirring by embryos and damping by gas 



drag is obtained as (jKobayashi et al.ll201 



rg I 
L0) 



e 2 = 



h 3 M (P VSt \ ow )TQ K 
2 7 /3 7r 6 7? 



(27) 



Using Eqs. (|T3|, pB}, and J27J) under the Stokes 
regime, we have the fragment mass mf at the low- 
mass end of collision cascade: 



where 



m f0 



mi = m i0 M 3/2 , 



221M lt &C L Q rj c / 3 



(28) 



1/3' 



.3/2 



For 



< h 



M ■ 



(29) 

Ida fc Nakazawa (1989) found 
that the dimensionless collision rate for e <C Hm 
is given by (P C ol,low) = 11.3V^S, where the coeffi- 



cient is determined by llnaba et al.l (|200lh . bince 
the atmosphere effectively enhances an embryo ra- 
dius for the accretion of bodies, the collision rate 



is modified to be (Inaba & Ikoma 2003) 



(-fcol) — (-PcoUow)l 



(30) 



We obtain the accretion rate of fragments by an 
embryo, M, from Eqs. and O as 



ifa 



F' 1 / 2 (F col i ow )C acc q 2 r?k 
mJ /18 (3M*) 2 / 3 



6/7 



(31) 



(32) 



Fragments with ni{ at the low-mass end of 
collision cascade that dominate the surface den- 
sity of solids E s are no longer disrupted by col- 
lisions and drift inward by gas drag. The drift 
velocity is given by 2rj 2 a/r and then the S s - 
reduction rate due to the radial drift is expressed 
as dTj s /dt — — 2(9/4 — (/)?7 2 E s /t with the assump- 
tion of E s cx a~ q . Since r of fragments with rrif is 
determined by Equations ([2"B)) and (|2"7)) . we have 



8 



(|Kobavashi et al.ll2010h 



dt 

Bf a — 



B fa E s M, (33) 

9 \ (-Pvs,low)^k?7Wk 



l- q 



2 4 / 3 37rM„C L 6QJ 



(34) 



Since fragments are later produced by embryo 
growth in an outer disk, the radial distribution 
depends on time in contrast to the assumption of 
S s oc a~ q . Nevertheless, the effect is negligible for 
embr yo growth unless the a tmosphere is consid- 
ered ( Kobavashi et al. 20101) . We discuss this ef- 
fect with the atmospheric enhancement in 2] and 

m 

We can now obtain the final embryo mass Mf a 
from M and S s for fragment accretion, similar to 
the case of planetesimal accretion. Integration of 
Equation (|31l) divided by Equation (TT51) results in 



M fa = 



/41A 



la 



36B fa 



42/41 



^36/41 
J s,0 ' 



(35) 



where we assume Mf a ^ Mo and S s .o 3> E s . For 
q = 3/2, we have 



M fa = 0.20 



a \ H7/164 



(5Au) 

/ice ^1 

30 gem -2 
Qh 



0.01 gem 



-3 



1/7 



36/41 



3.1 x 10 6 ergg- 1 



42/41 



M e . (36) 



Here, we adopted f lce and Si for the minimum 
mass solar nebula model. The weak dependence 
of Mf a on k implies that the overestimate of R c /R 
due to the power-law radial profile is insignificant, 
as discussed in Section [ 



For the case without an atmosphere . iKobavashi et al 
( 201dh derived a final mass for the fragment ac- 
cretion, 



M f = 0.14 



(— 

V5AU 



3/8 



Qh 



3.1 x 10 6 ergg 



/ice Si 

30 gem -2 

3/4 

M 



3/4 



(37) 



Eqs. (|36| and (1371) imply that the final masses 
increase due to the atmosphere, but the enhance- 
ment is insignificant; Mf a /Mf ~ 1.4-2 for 1- 
lOxMMSN. 



If we neglect the collisional enhancement due to 
atmosphere, the final mass M na is determined by 
the larger of M c and M f (Kobayashi et al. 2010). 
In the case with atmosphere, a final mass M a is 
also given by the larger of M ca and Mf a . The final 
mass M a is shown in Figs. [3H5] For the initial 
planetesimal radius tq = 10 km, M a is dominated 
by Mfa inside the point where the line of M a bends 
in Fig. [3] and by M ca outside. The final mass M a 
is determined only by Mf a for ro = 1 km (Fig. rj} 
and by M ca for ro = 100 km (Fig. [5]) in the range 
of interest. 



100 
10 
1 

0.1 

m 

10 

1 

0.1 

°M 

10 

1 

0.1 
0.01 




o 
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4 7 10 30 
Distance [AU] 

Fig. 3. — Embryo masses with (circles) and 
without (squares) atmosphere after 10 7 years for 
to = 4.2 x 10 18 g (ro = 10km), as a func- 
tion of distance form the central star. We set 
Si = 71gcm~ 2 (top), Si = 21gcm~ 2 (middle), 
and Si = 7.1gcm~ 2 (bottom). Solid lines in- 
dicate M a which is the larger of Af ca and Mf a 
for k = 0.01cm 2 g _1 . Dotted lines represent M na 
which is the larger of M c and Mf . Thin lines show 
M sn . 
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Fig. 4.— Same as Fig. [2 but for m = 4.2 x 10 15 g 
(radii of 1 km). 



0.1 
0.01 



4 7 10 30 
Distance [AU] 



Fig. 5.— Same as Fig. El but for m = 4.2 x 10 21 j 
(radii of 100 km). 



4. NUMERICAL SIMULATION 

Regarding the method of numerical sim ulation, 
we basicall y follow Kobavashi et al. ( 2010f ). The 
method of iKobavashi et al. ( 201dh is briefly ex- 
plained here. In the calculation, a disk is divided 
into concentric annuli and each annulus contains 
a set of mass batches. We set the mass ratio 
between the adjacent batches to 1.2, which can 
reproduce the collisional growth of bodies result- 
ing f rom TV-body simulatio n without fragmenta- 
tion ( Kobavashi et al] [2010) and the analytical so- 
lution of mass depletion due to collisional grinding 
(|Kobavashi fc Tana ka 2010). The mass and veloc- 
ity evolution of bodies and their radial transport 
are calculated as follows. 

- The mass distribution of bodies evolves 
through their mutual collisions that pro- 
duce mergers and fragments. The total mass 
of fragments ejected by a single collision is 



given by Equation @ and the remnant be- 
comes a merger. The collision rates between 
the bodies arc calculated from the formulae 



of llnaba et all (|200lh . 



The random velocities given by e and i of the 
bodies simultaneously evolve through their 
mutual gravitational interactions, gas drag, 
and collisiona l dam ping. The formulae of 
Ohtsuki et al] (|2002l) are applied to describe 
the changing rates of e and i. The gas- 
drag damping rates of e and i are described 
as functions of e, i, 77, and t according to 



Inaba et al.l (|200ll ). To determine r, we take 
into account Stokes and Epstein drag as well 
as a drag law with a quadratic dependence 
on velocity. For the collisional damping, 
both fragments and a merger resulting from 
a single collision have the velocity dispersion 
at the gravity center of colliding bodies. 
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In each annulus there is a loss and gain 
of bodies due to their inward drift. The 
number loss rate from an annulus is given 
by j (N (rnjVdiin / Aa)dm, where i>drift is the 
drift velocity of bodies, N(m)dm is the num- 
ber of bodies with mass ranging from m 
to to + dm in the annulus, and Aa is the 
width of the annulus. The bodies lost from 
each annulus are added to the next inner 
annulus. The drift velocity is given by 
(|Kobavashi et al.ll2010h 



Adrift 



2arj 



T 2 
' stop 



T 1 



f 2 
stop 



{2E + Kf 2 4 - 



9tt 2 



-e H 



(38) 

where E = 2.157, K = 1.211 and the dimen- 
sionless stopping time f stop = (e+i+rj) 
is adopted. 

In this paper, we add a collisional enhance- 
ment due to the atmosphere. Although the sim- 
ple power-law radial density profile of the atmo- 
sphere (Equation ©) is used for the derivation of 
final masses (M ca , Mf a ), the simulation incorpo- 
rates a more realistic profile pro vided by the for- 
mulae of llnaba fc Ikornal ((2003). The opacity of 
the embryo's atmosphere in their model is given 
by K = ftgas + /tgrj where K gas is the gas opacity, 
K gr is the opacity of grains having an interstellar 
size distribution, and / is the grain depletion fac- 
tor. Following Inaba & Ikoma, we adopt 

0.01 + 4/ cm 2 g _1 for T<170K, 
0.01 + 2/ cm 2 g" 1 for 170 K < T < 1700 K, 
0.01 cm 2 g" 1 for T > 1700 K. 

(39) 

The enhancement factor R e /R due to the atmo- 
sphere is shown in Fig. [2j 

We perform the simulations for embryo for- 
mation starting from a monodisperse mass pop- 
ulation of planetesimals of mass too and ra- 
dius ro with e = 2i — (2mo/M t ) 1 ^ 3 and p p 
1 g cm -3 around the central star of mass M Q 
with a set of eight concentric annuli at 3.2, 
4.5, 6.4, 9.0, 13, 18, 25, and 35 AU contain- 
ing E gas and E s for q = 3/2. To compute 
Qq, we use Equation (J5J) with Q 0s = 7.0 x 
H^ergg" 1 , /3 S = -0.45, Q n , = 2.1 erg cm 3 g~ 2 
p p = 1.19, and G,„ = 9 dBenz fe Asphaud Il99a 



IStewart fe Leinhardtl 12009^ We artificially ap- 
ply the gas surface density evolution in the form 



Egas — Eg a s, o exp(-£/T gaSidep ), where T gaSidep 
is the gas depletion timescale, which we set 
to r gaS! dop = 10 7 years. Assuming a constant 
E gas gives almost the same results for final em- 
bryo masses, because we consider time spans 

^ — -^gas,dep- 

Fig. [TJ] shows the embryo-mass evolution at 
6.4 AU for / = 0.01. Runaway growth initially oc- 
curs; embryo mass exponentially grows with time 
during the stage. The run away-growth timesc ale 
is proportional to ro/E s ,o ( Ormel et al. 2010al lb). 
j/gWhen the embryo masses exceed 0. 001-0. 01M®, 
oligarchic growth starts. Since massive embryos 
dynamically excite planetesimals, the reduction 
of planetesimals due to c ollisional fragmentatio n 
stalls the embryo growth (|Kobavashi et al. 201dl) . 
For E = 7.1 gem" 2 (MMSN), the fragmenta- 
tion limits the final mass to about Mars mass (~ 
0.1M ffi ) and the atmosphere is insignificant. Once 
embryo masses exceed the Mars mass, atmosphere 
substantially accelerates the embryo growth. For 
E > 21 gem" 2 (3 x MMSN), the atmosphere leads 
to further embryo growth. Nevertheless, embryos 
finally attain asymptotic masses. 

Results for these simulations are summarised in 
Fig.[3l where the embryo masses after 10 7 years are 
compared to analytical formulae for final embryo 
masses. Embryo masses finally reach M a inside 
5 AU (E = 7.1 gem- 2 ), 10 AU (E = 21 gem- 2 ), 
and 20 AU (Eo = 71 gem -2 ). However, embryos 
exceed M a inside 5AU for Eo = 71 gem -2 . This 
excess comes from the embryo growth through col- 
lisional accretion with bodies drifting from out- 
side, which effect we did not consider in the anal- 
ysis described in Section |3l To confirm the contri- 
bution from drifting bodies, we show the surface 
density evolution in Fig. [7l For Eo = 71 gem -2 , 
the surface density of solids increases after 2 x 
10 5 years. Since the drift timescale shortens in- 
ward, bodies from outside cannot raise the surface 
density unless embryos accrete them. Therefore, 
the increase in the surface density implies that em- 
bryo grows through the accretion of such bodies. 

The initial mass Too of planetesimals in the sim- 
ulations depends on their formation process, which 
is not well understood yet. We perform the em- 
bryo growth starting from different too (Figs. 2] 
and [5]). Small planetesimals are relatively easily 
fragmented due to low and quickly ground 
down to the low-mass end of collision cascade. 
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Fig. 6. — Evolution of embryo mass at 6.4 AU 
with m = 4.2 x 10 18 g (r = 10 km) for 
So = 71gcm~ 2 (lOxMMSN; top), 21 gem" 2 
(3xMMSN; middle), and 7.1 gem" 2 (MMSN; bot- 
tom). Solid lines show the case with atmosphere 
and dotted lines represent the result without at- 
mosphere. 

The resulting fragments with low e actively ac- 
crete onto embryos. For m Q = 4.2 x 10 15 g 
(t*o = 1 km) , embryos can reach a final mass M a 
in a relatively wide region inside 10 AU (MMSN), 
20 AU (3 x MMSN), and 30 AU ( 10 x MMSN). On 
the other hand, large initial plane tesimals de- 



lay th e runaway growth of embryos (jOrmel et al 



2010al lbh and the following oligarchic growth is 
also slower than that for small planetesimals be- 
cause embryos mainly accrete original planetes- 
imals rather than fragments with low e. For 
ro = 100 km, embryos attain the final masses 
only inside 4AU for 3 x MMSN and inside 6AU 
for 10 x MMSN, and embryos cannot reach final 
masses beyond 2.7 AU in the MMSN disk. In ad- 
dition, small bodies drifting from outside are ef- 
fectively captured by embryos and thereby em- 
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10 2 10 3 10 4 10 5 10 6 10 7 

Time [yr] 

Fig. 7. — The solid surface density evolution at 
3.2 AU. 

bryos exceed final masses M a inside 4AU for 
10 x MMSN. 

In the case without an atmosphere, initially 
larger planetesimals can form massive embryos. 
Since large planetesimals delay embryo growth, 
embryos made from lOOkm-sized initial planetesi- 
mals can reach 10 M® but the lo cation is only in- 



side 3 -4 AU even for 10 x MMSN (jKobavashi et al 
2010h . The case with the atmosphere shows a 



similar dependence of the final embryo masses 
on initial planetesimal mass. However, since the 
atmosphere accelerates embryo growth, embryos 
larger than 10 M e are produced inside 8-9 AU of 
a 10 x MMSN disk with lOOkm-sized initial plan- 
etesimals. 

While the final masses of embryos exceed 10M® 
for large initial planetesimals of ro > 100 km, em- 
bryos must reach the critical core mass within the 
disk lifetime T gaS! d p to form gas giant planets. 
The growth timescale is estimated to be M/M, 
where M is given by Equation (|16j) . The crit- 
ical distance a c inside which embryos can reach 
10M® is approximately obtained from the condi- 
tion M/M < T ga8idep with M = 10M e , 



9.6 



( T < 



dep 



20/39 



\ 10 7 years 

rp \ -V3 
100 km 



Si 



71 gem 



AU, 



23/39 



(40) 
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where we adopt m = 100 tuq and q 



3/2. 



For ro = 100 km, the massive disk with Si > 
70gcm~ 2 can form such large embryos around 
10 AU. In addition, we estimate a c ~ 5AU 
from Equation (gOJ) for a lOxMMSN disk with 
ro = 10 3 km. Indeed, the simulation with Si = 
71 gem -2 and ro = 100km shows embryos cannot 
reach 10M® beyond 5AU (see Fig-HJ. Therefore, 
the condition of lOxMMSN with r ~ 100km is 
necessary to form gas giants around 10 AU. 




7 10 30 
Distance [AU] 

Fig. 8.— Same as Fig. El but for Si = 71gcm- 2 
with m = 4.2 x 10 24 g (r = 1000km). The final 
mass M a with atmosphere is estimated to be larger 
than 200M®. 

We also give a constraint on /. For / < 0.01, a 
final mass is almost independent of / (see Fig. . 
This is because the gas opacity dominates over the 
grain opacity (see Equation (p9|V For / = 1, em- 
bryos at 3-4 AU become larger due to the capture 
of bodies drifting from outside, while final embryo 
masses in the outer disk are similar to the case 
without atmosphere. The condition of / < 0.01 
is therefore necessary for gas giant formation in 
the region 5-10 AU and such low / is acceptable; 
the depletion factor / should be much smaller 
than unity after planetesimal formation. In ad- 
dition, a low-op acity atmosph e re reduces the crit- 
ical core mass ( Mizuno Il98d: llkoma et all 12000 ; 
Hori fc Ikomall2010l) . 
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Fig. 9. — The final embryo masses for / = 0.0001 
(top), 0.01 (middle) and 1 (bottom), starting from 
lOxMMSN with m = 4.2 x 10 21 g (r = 100km). 
Lines and symbols are the same as in Fig. |3l but 
we apply k = 1cm 2 g _1 to derive M a for / = 1. 

5. DISCUSSION 

We derived final embryo masses analytically 
and numerically. They agree with each other quite 
well in the inner disk where the embryo forma- 
tion timescale is shorter than the nebula lifetime 
(~ 10 7 years). The analytical formula for final 
masses M a implies that initial planetesimal radii 
should be larger than about 3 x 10 3 km to form em- 
bryos with 10 M® at 5 AU in a MMSN disk. How- 
ever, the critical distance a c inside which embryos 
reach 10 M® within 10 7 years (Equation (|40[)) is 
estimated to be much smaller than 5 AU; a massive 
disk is likely to form gas giant planets. Embryos 
inside 5AU of a ~ 10 x MMSN disk exceed final 
embryo masses M a due to the accretion of small 
bodies drifting from outside. In spite of such fur- 
ther growth, embryos starting from small planetes- 
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imals cannot reach the critical core mass ~ 10 M®. 
In addition, further growth is insignificant beyond 
5 AU. The formulae for M a and a c suggest that 
initial planetesimals with ro ~ 50-700 km are nec- 
essary for embryos to reach 10 Mj, at 5 AU in the 
lOxMMSN disk. 

Inaba et al. (2003) performed similar simula- 
tions incorporating collisional fragmentation and 
enhancement due to the embryo's atmosphere 
and showed a planetary core with M > 10M® 
could be produced around 5 AU with tuq = 
4.2 x 10 18 g (r = 10km) for lOxMMSN. In 
our simulation, embryos cannot reach 10M® un- 
der this condition and larger planetesimals are 
necessary to form such massive embr yos beyond 
5AU. As Kobavashi fc Tanaka ( 2010l ) discussed, 
Williams fc Wetherilll ( 1994 ) underestimated the 
total ejecta mass produced by a single collision for 
cratering; Inaba et al. adopt ed the fragmentation 
model similar to theirs that IWetherill fc Steward 
(|l993l ) developed (see Fig. [1]). Erosive collisions 
shorten the depletion time of lOkm-sized plan- 
et eshniikj_Jj^^ by a factor of 4- 
5 (jKobavashi fc Tanakal 120101 ) and hence reduce 
final embryo masses. As seen from Eqs. 
and (|5fJ)) . final embryo masses M ca , Mf a in- 
crease with Q^; the results of Inaba et al. cor- 
respond to embryo masses for higher QJj,. Al- 
though we and Inaba et al. a pplied Q^> pro- 
vided by iBenz fc Asnhaud (|l999l) . porous bodies 
w ith r < 10km may have much lower Qn (e.g., 



sion may stop growth of smaller objects ([Okuzumi 
2009). A new scenario that allows one to over- 



2009; iMachii fc Nakamura 



Stewart fc Leinhardt 
20111 ). For initial planetesimals with radii > 
100 km, of slightly larger bodies determines 
final embryo masses and is almost entirely deter- 
mined by the gravitational binding energy; the 
uncertainty from their structure would be minor. 
Therefore, such large planetesimals are possible to 
produce cores for gas giant planets. 

The mechanisms of planetesimal formation are 
highly debated but, despite intensive effort, re- 
main fairly unknown. The formation through col- 
lisional coagulation in which dust smoothly grows 
to planetesimals with tq ~ 1 km face barriers: 
meter-sized objects should be lost to the central 
star as a result of g as drag (j Weidenschillingl 1 1 9 77 : 
Brauer et al J 120081). and further agglomeration of 
cm-sized objects upon collisio n is problematic be- 
cause of collisiona l bouncing ( Guttler et al. 2010t 
Zsom et al.l l201ol ). Moreover, the electric repul- 



come the barriers has been proposed recently: self- 
gravity of small particles accumulating in turbu- 
lent structures of gaseous disks fo rms large plan- 
etesimals of the order o f 100 km ( Johansen et al.l 
20071 : ICuzzi et al.l l2008h . Not only do such large 



planetesimals produce planetary cores exceeding 
the critical core mass to form gas giant planets, 
they may also be consistent with properties of mi- 
nor bodies in the solar system. Indeed, the initial 
planetesimals should be larger than 100 km to re- 
produce t he mass distribution of asteroids in the 
main belt (|Morbidelli et al.ll2009h . 

For large planetesimals, a final embryo mass 
given by Af ca is large enough to start core ac- 
cretion, while embryo growth is slow. If the ra- 
dial slope of surface density q = 3/2 like the 
MMSN model, a massive disk with lOxMMSN 
is necessary for embryos to reach the final mass 
around 10 AU. However, observations of proto- 
planetary disks infer their relatively flatter ra- 
dial distributions ove r several hundred AU (e.g., 
iKitamura et al.ll2002t ). In such a disk, dust grains 
accumulate in an inner disk due to radial drift dur- 
ing their growth, which i ncreases the solid s urface 
density in the inner disk ( Brauer et al.ll2008l) . The 
enhancement of solid surface density accelerates 
embryo growth and hence embryos may achieve 
the critical core mass in less massive disks. 

To form gas giants via core accretion, rapid 
gas accretion onto a core with ~ 10 Earth masses 
must occur prior to gas depletion. However, these 
cores migrate inward due to their exchange of an- 
gular momentum with the surrounding gas (Type 
I). From linear analysis, the characteristic orbital 
decay time of Earth-mass cores at several AU in 



the M MSN model is about 1 Myr (|Tanaka et al 
l2002h . Several processes to delay the timescale of 
Type I migration have been pointed out, for exam- 
ple, di sk surface density tra nsitions dMasset et al 



2006b), intrinsic turbulence ([Nelson fc Papaloizou 



2004), and hydrodynamic feedback (jMasset et al 



2006af ). There is still uncertainty about this es- 
timate of the migration time. Indeed, the dis- 
tribution consistent with observations of exoplan- 
ets can be reproduced only if the timescale of the 
type I migration is at least an order of magnitude 
longer than that derived from the linear analysis 
(llda fc Linl 2008). We should also investigate the 
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strength of such migration for the survival of cores 
of gas giant planets in our future work. 

6. SUMMARY 

In this paper, we investigate the growth of plan- 
etary embryos by taking into account, among oth- 
ers, two effects that are of major importance. One 
of them is collisional fragmentation of planetesi- 
mals, which is induced by their gravitational inter- 
action with planetary cores. Another effect is an 
enhancement of collisional cross section of a grow- 
ing embryo by a tenuous atmosphere of nebular 
gas, which becomes substantial when an embryo 
has reached about a Mars mass. 

The main results are summarized as follows. 

1. If the atmosphere is not taken into account, 
collisional fragmentation suppresses plane- 
tary embryo growth substantially. As a re- 
sult, embryos cannot reach the critical core 
mass of ~ 10M® needed to trigger rapid 
gas accretion to form gas giants. The fi- 
nal masses are about Mars m ass in a MMSN 
2010h . Embryo's at- 



disk ([Kobavashi et al 



mosphere accelerates the embryo growth and 
may increase the final embryo mass by up to 
a factor of ten. 

2. Planetary embryos attain their final masses 
asymptotically. We have derived the final 
mass analytically. The final mass of an em- 
bryo is predicted to be the larger of M ca and 
M fa , which are given by Eqs. jUJl and (|36) . 
respectively. These final masses are in good 
agreement with the results of statistical sim- 
ulations. 

3. Our solution indicates that an initial plan- 
etesimal radius ro>3x 10 3 km is necessary 
to form a planetary core with 10 M ffi at 5 AU 
in a MMSN disk. However, such initially 
large planetesimals delay embryo growth; a 
massive disk is required to produce massive 
cores within a disk lifetime. The analytical 
solution for the final mass and the embryo 
formation time show that planetesimals with 
an initial radius of tq ~ 50-700 km are likely 
to produce such a large planetary core within 
a disk lifetime at 5 AU for 10 x MMSN. 



4. The embryo growth depends on the disk 
mass, initial planetesimal sizes, and the 
opacity of atmosphere. We have performed 
statistical simulations to calculate the final 
embryo masses over a broad range of param- 
eters. We took the surface density of solids 
at 1 AU in the range of Ei = 7.1-71gcm -2 
(1-10 x MMSN), initial planetesimal radius 
ro = 1-1000 km, and the grain depletion 
factor / in planetary atmosphere between 
/ = 10~ 4 -1. We found that planetary em- 
bryos can exceed IOMq within 8-9 AU for 
lOxMMSN, r = 100km, and / < 0.01. 
Other sets of parameters cannot produce 
massive cores at 5-10 AU. For example, em- 
bryo's mass can reach 6M ffi for ro = 10 km 
only inside 4AU. Therefore, we conclude 
that a massive disk (~ 10 x MMSN) with 
ro ~ 100 km and / < 0.01 is necessary to 
form gas giant planets around 5-10 AU. This 
condition for large embryo formation is in- 
dependent of the material strength and/or 
structure of bodies, because QJj of 100km- 
sized or larger bodies is largely determined 
by their self-gravity. 



We thank Chris Ormel for helpful discussions 
and the reviewer, John Chambers, for useful com- 
ments on the manuscript. 

REFERENCES 

Adachi, I., Hayashi, O, & Nakazawa, K. 1976, 
Progress of Theoretical Physics, 56, 1756 

Benz, W., & Asphaug, E. 1999, Icarus, 142, 5 

Bodenheimer, P., & Pollack, J. B. 1986, Icarus, 
67, 391 

Bottke, W. F., Durda, D. D., Nesvorny, D., 
Jcdicke, R., Morbidelli, A., Vokrouhlicky, D., 
& Levison, H. F. 2005, Icarus, 179, 63 

Brauer, F., Henning, T., & Dullemond, C. P. 2008, 
A&A, 487, LI 



Chambers, J. E. 2006, ApJ, 652, L133 



15 



Chambers, J. 2008, Icarus, 198, 256 

Cuzzi, J. N., Hogan, R. C, & SharifT, K. 2008, 
ApJ, 687, 1432 

Dohnanyi, J. S. 1969, J. Geophys. Res., 74, 2531 

Fujiwara, A., Kamimoto, G., & Tsukamoto, A. 
1977, Icarus, 31, 277 

Greenzweig, Y., & Lissauer, J. J. 1992, Icarus, 100, 
440 

Guttler, C, Blum, J., Zsom, A., Ormel, C. W., & 
Dullcmond, C. P. 2010, A&A, 513, A56 

Hayashi, C. 1981, Progress of Theoretical Physics 
Supplement, 70, 35 

Holsapplc, K. A. 1993, Annual Review of Earth 
and Planetary Sciences, 21, 333 

Hori, Y., & Ikoma, M. 2010, ApJ, 714, 1343 

Housen, K. R., Schmidt, R. M., & Holsapple, K. A. 
1991, Icarus, 94, 180 

Ida, S., & Lin, D. N. C. 2008, ApJ, 673, 487 

Ida, S., & Makino, J. 1993, Icarus, 106, 210 

Ida, S., & Nakazawa, K. 1989, A&A, 224, 303 

Ikoma, M., Nakazawa, K., & Emori, H. 2000, ApJ, 
537, 1013 

Inaba, S., & Ikoma, M. 2003, A&A, 410, 711 

Inaba, S., Wetherill, G. W., & Ikoma, M. 2003, 
Icarus, 166, 46 

Inaba, S., Tanaka, H., Nakazawa, K., Wetherill, 
G. W., & Kokubo, E. 2001, Icarus, 149, 235 

Inaba, S., Tanaka, H., Ohtsuki, K., & Nakazawa, 
K. 1999, Earth, Planets, and Space, 51, 205 

Johansen, A., Oishi, J. S., Mac Low, M.-M., 
Klahr, H., Henning, T., & Youdin, A. 2007, Na- 
ture, 448, 1022 

Kenyon, S. J., & Bromley, B. C. 2004, AJ, 127, 
513 

Kenyon, S. J., & Bromley, B. C. 2008, ApJS, 179, 
451 



Kenyon, S. J., & Bromley, B. C. 2009, ApJ, 690, 
L140 

Kitamura, Y., Momosc, M., Yokogawa, S., 
Kawabe, R., Tamura, M., & Ida, S. 2002, ApJ, 
581, 357 

Kokubo, E., & Ida, S. 1996, Icarus, 123, 180 

Kokubo, E., & Ida, S. 1998, Icarus, 131, 171 

Kokubo, E., & Ida, S. 2000, Icarus, 143, 15 

Kokubo, E., & Ida, S. 2002, ApJ, 581, 666 

Kobayashi, FL, & Tanaka, H. 2010, Icarus, 206, 
735 

Kobayashi, H., Tanaka, H., Krivov, A. V., & In- 
aba, S. 2010, Icarus, 209, 836 

Levison, H. F., Thommes, E., & Duncan, M. J. 
2010, AJ, 139, 1297 

Machii, N., & Nakamura, A. M. 2011, Icarus, 211, 
885 

Masset, F. S., D'Angelo, G., & Kley, W. 2006, 
ApJ, 652, 730 

Masset, F. S., Morbidelli, A., Crida, A., & Fer- 
reira, J. 2006, ApJ, 642, 478 

Mizuno, H. 1980, Progress of Theoretical Physics, 
64, 544 

Morbidelli, A., Bottke, W. F., Nesvorny, D., & 
Levison, H. F. 2009, Icarus, 204, 558 

Nelson, R. P., & Papaloizou, J. C. B. 2004, MN- 
RAS, 350, 849 

Ohtsuki, K. 1999, Icarus, 137, 152 

Ohtsuki, K., Stewart, G. R., & Ida, S. 2002, 
Icarus, 155, 436 

Okuzumi, S. 2009, ApJ, 698, 1122 

Ormel, C. W., Dullcmond, C. P., & Spaans, M. 
2010a, ApJ, 714, L103 

Ormel, C. W., Dullcmond, C. P., & Spaans, M. 
2010b, Icarus, 210, 507 

Stevenson, D. J. 1982, Planet. Space Sci., 30, 755 

Stewart, G. R., & Ida, S. 2000, Icarus, 143, 28 



16 



Stewart, S. T., & Leinhardt, Z. M. 2009, ApJ, 691, 
L133 

Takagi, Y., Mizutani, H., & Kawakami, S.-I. 1984, 
Icarus, 59, 462 

Tanaka, H., Takeuchi, T., & Ward, W. R. 2002, 
ApJ, 565, 1257 

Tanigawa, T., & Ohtsuki, K. 2010, Icarus, 205, 
658 

Thommcs, E. W., Duncan, M. J., & Levison, H. F. 
2003, Icarus, 161, 431 

Weidenschilling, S. J. 1977, MNRAS, 180, 57 

Wcidenschilling, S. J. 2005, Space Sci. Rev., 116, 
53 

Weidenschilling, S. J. 2008, Physica Scripta Vol- 
ume T, 130, 014021 

Weidenschilling, S. J., Spaute, D., Davis, D. R., 
Marzari, F., & Ohtsuki, K. 1997, Icarus, 128, 
429 

Wetherill, G. W., & Stewart, G. R. 1989, Icarus, 
77, 330 

Wetherill, G. W., & Stewart, G. R. 1993, Icarus, 

106, 190 

Williams, D. R., & Wetherill, G. W. 1994, Icarus, 

107, 117 

Zsom, A., Ormel, C. W., Guttler, C, Blum, J., & 
Dullemond, C. P. 2010, A&A, 513, A57 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



17 



